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We present the status of a lattice calculation for the K — > tttv matrix elements of the AS — 1 effective weak 
Hamiltonian, directly with two pion in the final state. We study the energy shift of two pion in a finite volume both 
in the /=0 and 1—2 channels. We explain a method to avoid the Goldstone pole contamination in the computation of 
renormalization constants for AJ=3/2 operators. Finally we show some preliminary results for the matrix elements 
of AI—1/2 operators. Our quenched simulation is done at f3—Q.O, with Wilson fermions, on a 24^ x 64 lattice. 



1. General Strategy 

Kaon weak decay amplitudes can be described in 
terms of matrix elements (ME's) of the AS'^1 ef- 
fective weak Hamiltonian. H^^^^ is written as a 
linear combination of a complete basis of renormal- 
ized local operators (OP's) Qi{fi), where fj, is the 
renormalization scale. The most relevant contribu- 
tions in the computation of the weak amplitudes 
■^1=0,2 and of e'/e are given by the K^tttt ME's of 
Q~, Qe, Qt and Qg (see below and Ref. [|l|,||). 
In order to compute these ME's from lattice QCD, 
one has to renormalize the bare (divergent) lattice 
OP's. Two methods are possible on the lattice: 
1) ComputejO|Qi(/Lt)|iC) and {Tr\Qi{^j,)\K) and then 
derive {TnT\Qi{iJ,)\K)i=Q^2 using soft pion theorems [ 
^. In this case, besides the difficulties in controlling 
chiral behaviour, a major problem is the inclusion 
of higher order contributions in the chiral expansion 
which should reconstruct the final state interaction 
(FSI). In fact, in this approach, only K ^ tttt ME's 
at lowest order in ChPT can be obtained (see Refs. [ 

1) . „ 

2) Compute directly {TTn\Qi{^)\K) i=o,2 - The main 
difficulty in this case is the relation between ME's in 
a (Euclidean) finite volume and the corresponding 
physical infinite volume ones. Even if, in principle, 
it is possible to take into account FSI exactly [ ^, |^ , 
in practice these methods are numerically very de- 
manding. So, with present computing power, only 
simulation at unphysical kinematics [ |^ may be 
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achieved, and ChPT is still needed to extrapolate 
to the physical point. 

In this work we choose the second strategy. To 
extract ME's we place the local OP Qi in the origin 
and we use three local interpolating fields (e.g. the 
pseudoscalar density P"*" — ^p'^j^ip^), one at time 
t2 — 10 which annihilates one tt in the two pion 
state and one at time tx — 54 = —10 which creates 
a K (both with momentum zero). As explained 
below t2 {tx) must be chosen large enough in mo- 
dulus for the two pion state (the kaon state) to be 
already asymptotic. The third one annihilates at 
time ti (not fixed) the second tt with momentum 
either or 27r/L. As shown in [ ||], in the limit 
r/2>ti>t2>0, T>tK>r/2 we have 

{0\T[7:^{h)%2{t2)Q^{0)k{tK)]\0)v ^^^^ _^ 

GTr{ti)GTr{t2)GK{tK) ^ 

\{t:7:\Q^\K)\ COS S{W) e''"'"'''''"'+ oQ) 

where G^{t) = (0|T[7r(t)7rt(0)]|0)y, = \{0\tt\tt)\ 
(the definition of Gxit) and Zk is analougous), 
Itttt) is the lower two pion state with the same mo- 
mentum of ni{ti). AE — W — 2Mt^ is the energy 
shift of two pions in a finite volume (FV) while 
6{W) is the strong interaction phase (which de- 
pends on the isospin and on the energy of the two 
pion state). In order to extract the physical ME, 
one needs to compute FV corrections, which are 
represented by the correction to the ME (indicated 
above by the term 0{1/L)) and by the FV energy 
shift. We focus now on this last issue. Theoretical 
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Figure 1: a. G4.(t)/G.(t)" 

b. effective mass of G2Trs{t)/G,r{t)^ 

predictions [0, H) show that the factor exp{— AEt2) 
should give relevant corrections: of order 5 10% in 
the A/ — 3/2 channel and of order 25% (or larger) 
in the A/ = 1/2 one. 

2. Finite Volume Energy Shift 

To extract AE we study the four point function 
Ginit) = {0\T[TT4,{t)TT3{t)TTl{0)TTliO)]\0)v. Siuce the 
computation of the / = correlator require the 
computation of T Si propagators, we use 6*4^ (<) 
only for 1 = 2. We have 7r2(0)^i(0) = P(0)P(0) 
while Tr^{t)n3{t) = X^x *5xp iqx^^y -P(y, 
with q = 0, 27r/i and with the appropriate flavour 
structure. Due to the possibility of propagation 
both of a two pion state in one temporal direction 
and of a state composed by two single pion states 
in opposite directions, in the region T ^ t ^ one 
has 

C^^ cosh(W^(t - f )) + 
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where C,r7i- and C,r depend on the ME's of the 
interpolating fields and on the energy shift. We 
thus make a fit with Aexp{—AEt) in the region 
r/2 > i > (see Fig. 0.a). 

In the 7 = channel we use instead G2-Ks{t) = 
{0\T[S(t)Tfl{0)TTl{0)]\0)v where S is the scalar den- 
sity summed over the space and 7r2(0)7ri(0) = 
n(0) Ex Pi'^^ 0)' with n equal either to P or to Aq 
(the temporal component of the axial current). We 
have computed only the connected part, because we 
expect the disconnected one to be negligible since 
it is suppressed by the OZI rule. To check this as- 
sumption, a computation of the disconnected part is 
presently under way. We observe empirically that, 
for IL = Ao, the noise is smaller and the plateaux 
start at smaller times (already t ~ 10 seems suf- 
ficients, see Fig. |l|.b). Results for both 7 = 2 and 
7 = channel are shown in Table |^. We compare 



Table 1: Results in lattice units multiplied by 10"^ 



numerical results both with Liischer formula 
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(where ag~°'^ is the tree level result in ChPT) and 
with its version AEbg in 1-loop quenched ChPT [ 
H which, in the 7 = case, shows enhanced finite 
volume correction. For 7 = 2 there is very good 
agreement between numerical values and theoreti- 
cal predictions. For 7=0 there is a striking disagree- 
ment for the two heaviest masses while in the light- 
est case, also due to the large errors, they are com- 
patible. Moreover the dependence with the mass 
of the pion is different to the one predicted. This 
could be due to the presence, for large unphysi- 
cal quark masses, of a stable scalar particle below 
the two pion state. In this case the behaviour of 
the shift measured should have only exponentially 
small correction in the volume. Instead the "gen- 
uine" FV energy shift of the two pion state has a 
power dependence on the volume. To clarify the sit- 
uation we are thus currently varying the volume of 
our simulations and trying to reach smaller masses. 

3. Renormalization and Goldstone Pole 

The renormalization of A7 = 3/2 4-fermion 
OP's, which mix only with OP's of the same di- 
mension, may be achieved with the method of 
ref. [ ^. This requires the existence of a win- 
dow Aqcd <C /X <C 1/a, where the first inequal- 
ity serves to avoid pure non-perturbative effects 
like the coupling to the Goldstone boson. In 
some cases [ |lO|, this makes the chiral limit 
ill-defined and causes a systematic error in the 
extrapolation. It is easy to show that for par- 
ity conservig OP's this coupling gives both a sin- 
gle and a double pole, while in the parity violat- 
ing (PV) case only a single pole is present. For 
the sake of illustration, we explain the procedure 
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Figure 2: a. Da-iiMlp') b. (Z(M)Z(Mo)i7(^i, ^o))33 

to remove this unwanted contribution only in the 
PV case. Consider the amputated Fourier trans- 
form (at equal external momenta p) of the Green 
functions Gi = {^p{xi)^l}{x2)Oi{0)^p{x3)^p{x4)) . By 
projecting this array onto the possible indipendent 
Dirac structure we obtain a 5 x 5 matrix D{M^,p^) 
from which the matrix of RC's is obtained [ |l^ by 



j2, where Zj, is the quark 



field RC. The pole is present e.g. in -D33, which 
corresponds to the OP O3 = SP - PS. This is 
evident from Fig. |^.a. To subtract it, it is possi- 
ble either to fit each element of D directly with 
+ B{p'^)/Ml + C{p^)M^ or to construct the 
combination 



D{Ml,p^)Ml ~ D{Mlp^)Mi 



(1) 



(for non-equal masses), which automatically cancels 
the pole [ and fit it with A{p^) + C{p^){Ml + 
M|). Finally, by using either A{p'^) or A{p'^), we ob- 
tain two determination of the RC's perfectly com- 
patible with each other. In Fig. |[b the behaviour 
of {Z{ij,)ZIij,q)U{/j, no))^^ (where [/(/i,/io) is the 
renormalization group evolution between fio and fi 
at NLO) is shown with fio = 1.5 GeV, before and 
after the subtraction. It is clear that there is a win- 
dow for 2.4 GeV > fj. > 1.4 GeV where Z follows 
the expected renormalization group evolution. 

4. Very briefly on A/ = 1/2 ME's 

AI = 1/2 OP's mix, through penguin con- 
tractions, also with lower dimensional OP's with 
power divergent coefficients. In this case a non- 
perturbative subtraction is needed (in both strate- 
gies: K TTTT or K ir). Here we will present 
very preliminar results for with the purpose of 
showing that, for the first time, a signal has been 
observed (further details can be found in Ref. [ Q). 
Q~ enters, together with , in the computation 
of Re^7=o ■ They are defined as 
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Figure 3: (ttttIQ \K) bare on 340 configurations 

=s'^f:uuYld±sj^duj'lu-{u^c). (2) 

In our simulation the charm quark is propagat- 
ing. Due to the GIM mechanism, the subtraction is 
thus implicit in the difference of penguin diagrams 
with an up quark and a charm quark inside the 
loop. In Fig. ^the penguin contractions of the bare 
0~ are shown, in a kinematical configuration in 
which Mkh = 2MT^a = 0.72. It is interesting, even 
though only indicative, to note that the ratio of 
bare OP's {TnT\Q-\K)i=o/ {'ktt\Q+\K)i=2 « 9. In 
fact, although affected by a huge statistical error 
and incomplete for the lacking of the RC's and of 
part of the contractions, this result shows that pen- 
guin contractions are of the right order of magni- 
tude needed to explain the A/=l/2 rule (remem- 
ber that the ratio of Wilson coefficients |C~(/i = 
2GeV)/C+(/i = 2GeV)| « 2). Concerning 0+, 
penguin contractions gives a value close to zero (but 
again with large errors). In order to reduce the sta- 
tistical errors, a study to find the best source for 
the two pion state is presently under way. 
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